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Abstract 

We propose, in a heuristic way, a relativistic modified gravity model as an alternative to particle 
dark matter at galactic scales. The model is based on a postulated preferred time foliation described 
by a dynamical scalar field called the "Khronon". In coordinates adapted to the foliation it 

^^ , appears as a modification of general relativity violating local Lorentz invariance in a regime of 

weak gravitational fields. The model is a particular case of non-canonical Einstein-aether theory, 
but in which the «ther vector field is hypersurface orthogonal. We show that this model recovers 
the phenomenology of the modified Newtonian dynamics (MOND) in the non-relativistic limit, 

^sq ■ and predicts the same gravitational lensing as general relativity but with a modified Poisson-type 

potential. 

a^, PACS numbers: 95.35.+d,95.36.+x,04.50.Kd 
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I. INTRODUCTION 



Modifying gravity in the regime of weak gravitational fields or accelerations, is the leit- 
motiv of MOND — the Modified Newtonian Dynamics [jj-y]. The MOND paradigm is very 
successful at solving the problem of the fiat rotation curves of galaxies without the need of 
dark matter (see [^ for a review), and at explaining the correlations between the distribu- 
tion and dynamics of ordinary matter versus dark matter |5|, |6|. However, it is based on 
a non-relativistic formula and needs to be generalized in order to be applied in cosmology 
and to the gravitational lensing observed at galaxy cluster scale |7H9|, which are part of the 
great successes of the standard cosmological paradigm A-CDM based on particle cold dark 
matter [lO| and the cosmological constant A. 

We adopt for MOND the modification of the non-relativistic Poisson equation [llj. 
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where p is the density of ordinary (essentially baryonic) matter, is the gravitational po- 
tential, g = — V0 is the gravitational field and g = \g\ its ordinary Euclidean norm. The 
MOND function /i depends on the single argument g/ao, where oq = 1.2 x 10~^°m/s^ is a 
constant acceleration scale. In the limit of weak gravitational fields i.e. g '^ oq, the MOND 
function is linear, p{g/ao) = g/ao, while in the strong field regime g ^ a^ (though non 
relativistic), p tends to one so that (II. ip reduces to the usual Poisson equation. 



Several relativistic extensions of MOND have been proposed fse e |12 



^v.,..xi^x x...i^^x,.o^xv. wv^v...oxv^..o ^^ ^,^^^,^ ..c*,,.. >^^.... ^^^^^^~^^ , ^~^ y i ...^,. | 13[| for reviews). 
The Tensor- Vector-Scalar (TeVeS) theory of Bekenstein and Sanders |14J-|l6| extends general 
relativity with a time-like vector field and one scalar field. Einstein-aether theories IT], llSl 
involve a unit time-like vector field which is non- minimally coupled to the metric, and were 
originally used to describe the phenomenology of local Lorentz invariance violation. Those 
theories, when modified to involve non-canonical kinetic terms, have been shown to provide 

The vector field postulated 

The 



interesting examples of relativistic MOND theories [19|, |20 



in TeVeS is in fact analogous to the vector field in Einstein-aether theories |21l. 122 



cosmology in TeVeS and Einstein-^ther theories has been extensively studied 23|-|26| . The 
parametrized post-Newtonian parameters, included the preferred frame parameters crucial 
in the presence of a vector field, have been worked out for both the canonical Einstein-^ther 
theory J27| and TeVeS | 28l] . More recent proposals for relativistic MOND include a bimetric 
theory of gravity 29|, |30|, a refinement of TeVeS u sing a Galileon "fc-moufiage" to prevent 
deviations from general relativity at small distances |3l| , and a non-local extension of general 
relativity |32|. One can also invoke some new exotic properties of the dark matter rather 



than modifying gravity |33l . |34 



On the other hand, much interest focused recently on the possible violation of the local 
Lorentz invariance in a completion of general relativity at high energy. This Hof ava-Lifshitz 
approach |35| is motivated by the construction of a power-counting renormalizable quantum 
theory of gravity. The Lorentz invariance violation results from a preferred time foliation 
which is described in a 4-dimensional covariant formulation by a dynamical scalar field 
sometimes called the "Khronon" |36|. A consistent extension of the (non-projectable version 



of the) Hofava gravity has been proposed |36l . |37| to avoid some problems regarding the 
stability of the theory. In this extension the added terms crucially involve the acceleration 
of the congruence normal to the preferred family of hypersurfaces. It has been shown that 
the low energy limit of the extended Hofava gravity is equivalent to an Einstein-^ther theory 
in which the vector field is hypersurface-orthogonal J36l . l38l . 



In the present paper, we revisit some works on Einstein-aether theories with non-canonical 
kinetic terms (notably [19|, l20|), but in the framework of a preferred time foliation, inspired 



by these recent approaches to quantum gravity. We propose a relativistic MOND theory 
which modifies general relativity by terms depending on the acceleration in a way somewhat 



similar to the extended Hofava gravity |35l-l37l|. However, our motivation is quite different 



from that in |35|-|37]. as it is purely classical, and we are not concerned with the problem 
of quantizing gravity. In particular, the violation of Lorentz invariance will be effective at 
low energy, for very weak gravitational fields. Our aim is to provide a phenomenologically 
viable alternative to the particle cold dark matter (at least at galactic scales). 

The modification we propose consists of introducing in the Lagrangian of general relativity 
a single free function of the norm a of the acceleration of the congruence normal to the 



preferred hypersurfaces. (Unlike in [36|, |37|, we are not considering a finite number of terms 
including this acceleration, each one with a coupling constant scaling with the Planck mass, 
but a fully non-linear function.) The theory will be a particular case of non-canonical 
Einstein-aether theory but in which the aether field is hypersurface orthogonal. We shall 
present the theory in both a 3+1 decomposition and 4-dimensional covariant form, and 
shall explicitly prove the equivalence of the two formulations. In agreement with general 



studies on non-canonical Einstein- aether theories |19l . |20| , we will show that it is possible, 
by specifying the asymptotic behaviour of the function of the acceleration when a <^ oq, to 
recover the Bekenstein-Milgrom equation (II. ip in the non- relativistic limit. The theory will 
also be viable regarding the gravitational lensing of photons, which will be the same as in 
general relativity but with a Newtonian potential obeying the modified Poisson equation]^ 
The plan of this paper is as follows. In Section [TTl we present the 4-dimensional version of 
the theory, introducing the dynamical Khronon field defining the space-time foliation, and 
investigating the resulting field equations. In Section IIIII we switch to an alternative 3+1 
point of view, for which the time coordinate coincides with the Khronon field, and where 
the new terms associated with the acceleration in the action may be interpreted as Lorentz 
invariance breaking terms due to a preferred- frame effect. In Section [TVt we investigate the 
non-relativistic limit and the conditions under which one can retrieve the non-relativistic 
MOND equation (II. ip . Section |V] contains some concluding remarks. Finally, Appendix \K\ 
is devoted to the proof of equivalence between the Ad covariant formulation and the 3+1 
formalism, while Appendix |B]presents the canonical decomposition of the energy-momentum 
tensor associated with the Khronon field when interpreted as a matter field. 



II. COVARIANT FORMULATION 

In the 4-dimensional covariant formulation, we extend general relativity by means of the 
additional Khronon scalar field denoted by r. We define a preferred foliation of space-time 
according to the constant-r hypersurfaces, so we will require r to be smooth and free of 
extrema. In Section IIIII below, we consider a frame adapted to this foliation in which the 
time coordinate t will coincide with r, but in the present Section we investigate the fully 
covariant formalism. 

As in the usual 3+1 formalism, we introduce the unit timelike vector n^, oriented to- 
ward the future (n" > 0), orthogonal to the constant-r hypersurfaces, and normalized to 



^ While this work was in preparation, a paper based on similar ideas appeared on the archives 39 1. However, 
the field equations in that paper are in disagreement with ours. 



g^uTi^^n^ = — lo Thus we have n^ = —Nd^r, where A^ will later correspond to the ordinary 
lapse in adapted 3+1 coordinates and reads 

N= ^ ^ ^=. (2.1) 

The vector n^ is invariant by reparametrization r -^ f^r) where f (r) denotes any strictly 
increasing function. We also define the projector on the constant-r hypersurfaces by 7^j, = 
Ofiiy + n^Uy, and as usual the associated projected covariant derivative D^ (for instance 
D^ = 'Jn'^diy when acting on scalars; see Appendix |A] for more notation). The acceleration 
of the congruence orthogonal to the constant-r hypersurfaces is given by a'^ = n'^'Vyn'^ or, 
equivalently, in terms of the projected spatial derivative operator, 

a^ = D^\nN. (2.2) 

We will see in Section HVl that, in the non-relativistic limit and in adapted coordinates, the 
spatial components of the acceleration reduce to minus the ordinary gravitational accelera- 
tion. Hence observers which would follow the congruence in their motion with 4- velocity n'^ 
are essentially un-accelerated observers in coordinates adapted to the constant-r hypersur- 
faces. Thereafter we will call such observers the fiducial observers. 

In the present paper, motivated by the construction of a relativistic MOND theory mod- 
ifying gravity in a regime of small accelerations, we investigate a specific extension of the 
Einstein-Hilbert action for gravity by a function of the norm of the (space-like) acceleration 
(jMD, defined by 

a = ^lf,u a'^a'' . (2.3) 

In keeping with MOND we shall assume that such an extension becomes effective only 
when gravity is weaker than the MOND scale Oq (see Section HV]) . Well above Oq, but 
still remaining purely classical, we will assume that the standard general relativity with a 
cosmological constant is recovered. 

We consider a Lagrangian density C[g^^, r, \1/] depending on three dynamical fields, the 
metric g^u, the Khronon field r associated with the time foliation, and the matter fields 
collectively denoted by ^. Specifically we choose 



R-2f{a) +C4g^,,<il], (2.4) 
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where R is the curvature scalar, C^ denotes the Lagrangian of the matter fields, and /(a) 
is our postulated modification, in the form of an unspecified function of the norm of the 
acceleration (12. 3p . Only the asymptotic behaviour of /(a) when a — )■ oo and when a — )■ 
will be specified in Section HVl The factor in front of /(a) is chosen so that a constant term 
in /(a) will correspond to the usual definition of the cosmological constant. 

The specific choice for the Lagrangian (12. 4 p corresponds to one particular term inves- 
tigated in MOND-motivated non-canonical versions of Einstein-aether theories: see for in- 
stance the C4-term in the equation (30) of [20]. This term, among the simplest kinetic terms 
considered in Einstein-aether theories, is sufficient for our purpose of reproducing the MOND 
dynamics. It is also the only one to cancel out in an homogeneous and isotropic Universe 



^ We adopt the sign conventions of [40], and notably the ( — h ++) signature for the metric. Greek indices 
are space-time indices, while Latin are space indices. We pose G = c = I through most of the paper. 



(see e.g. [26|]). Note that the expression (12. 2 p we employ here for the acceleration is valid 
only in the hypersurface-orthogonal case, and implies for instance that Dfj,av = D^Qf^. 

The function /(a) is the only term depending on the Khronon field r in the action (12.41) . 
In adapted coordinates, this term will be interpreted as a Lorentz invariance breaking term 
associated with the preferred frame effect. Here we assume that the matter Lagrangian £j^ 
is standard, i.e. given by the standard universal coupling to the metric g^^, without depen- 
dence on the r-field. Apart from that, we do not make any restriction on the Lagrangian of 
the matter fields. 

Varying the action with respect to the metric yields a modified Einstein field equation, 

G^'" + f{a)9^'' + 2n^n^Vp [x{a)af] - 2x(a)a^a'^ = Stt T"'' , (2.5) 

whose trace gives 

R - 4/(a) + 2V/, [x(a)a^] + 2x{a)a'^ = -Svr T . (2.6) 

Here G'^'^ = R^'^ — ^gf^'^R is the Einstein tensor, T^^'^ = (^j ^—g^bC^jbg^y is the stress-energy 
tensor of the matter fields, and we have introduced the convenient notation 



X(a) 



2a 



(2.7) 



with /'(a) = d//da. Such notation is motivated by the non-relativist ic MOND limit, in 
which x(a) will appear as a gravitational analogue of the electric susceptibility of the (phan- 



tom) dark matter medium in the sense of [41 



Varying next with respect to the scalar field r leads to an independent equation, called 
the r-equation in the following, and given by (see J38[)o 



Nin^'Vyixa^) - xaV - x {K"" - -/""K) a. 







(2.8) 



where the extrinsic curvature tensor, which is symmetric, is defined by K^"^ = V^n^ + n^a" 
(with trace given by i^ = V^ra^). An alternative, simpler form of this equation, reads as 



n^V.(xa^) 



+ XCi^a''K^ 



/lU } 



(2.9) 



where we denote / = n'^Viyf. Note that this r-equation is of fourth order in derivatives of 
the Khronon field r. However, as will become clear in Section IIIIl if we select a coordinate 
system for which t = t, the r-equation will be equivalent to an equation of first order only 
in time derivatives of geometrical quantities. 

When we rephrase the model by interpreting r as a matter field, i.e. rewriting the 
modified Einstein equation in the form G'^'^ = 87r(T^^ + T^'^), where the stress-energy tensor 
of this matter field (say a "dark energy" fluid) reads 






fg^-" + 2n^n''Vpixan - 2xa^'a 



(2.10) 



■^ In the case where the aether vector is not hypersurface orthogonal, the quantity defined by K^'^ = V^n'^ + 
n^ay is not symmetric anymore (and does not warrant to be called an extrinsic curvature). Varying in this 



case the action with respect to n^ , one obtains instead of the scalar equation 



the vectorial equation 



n^V^ixa^ 



2 u 

Xa n 



x[K'"'-l'"'K)a,^Q, 



which has three independent components. Varying with respect to the metric, one also gets a different 
modified Einstein equation (|2.5p . with one additional term in the left-hand side given by —ixapU^^A"'^ 
where we pose Af"" = Rf"" - K"''. 



the r-equation is equivalent to the conservation of this associated energy-momentum tensor, 
VuT!^'^ = 0. Thus, given the modified Einstein equation (12.51) and the contracted Bianchi 
identity Vi^G'^" = 0, we see that the r-equation is in fact equivalent to the conservation 
of the matter energy-momentum tensor V^T'^" = 0, which is itself the consequence of the 
matter field equations and the scalarity of the action (see |38|] for discussion). In Appendix 
[B| we will give the canonical decomposition of the stress-energy tensor (12.101) . 

III. COORDINATES ADAPTED TO THE FOLIATION 

In this Section, we express the model in an equivalent 3+1 formalism, choosing a coordi- 
nate system adapted to the foliation. Note that this re-writing is not possible in a generic 
Einstein- aether theory like in [l9|, |20|. The time coordinate will coincide with the Khronon 



field: t = T and thus n^ = (— A^, 0) with A^ = l/y— (7°". Introducing the shift vector 
Aj = Qoi and the spatial metric 'jij = Qij, we have the usual 3+1 form for the metric, A^ 
playing the role of the lapse: 

ds^ = -{N^ - N,N')dt'^ + 2Nidtdx' + -fijdx'dx^ . (3.1) 

In particular n^ = ;^(1, — A^*). All spatial (Latin) indices are raised and lowered using the 
spatial metric 7jj and its inverse 7*-' = g^^ + N^N^ /N"^ (such that Y'^'^jk = ^l) which coincides 
with the spatial components of the 4-dimensional projector 'f^'^. 

The 4-dimensional Lagrangian (12. 4p now becomes, after discarding an irrelevant total 
divergence term (with y/—g = N^) , 
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n + K,,K'^ -K'- 2/(a) + C^[N, N„ 7,,, ^] . (3.2) 



Here 7^ denotes the 3-dimensional scalar curvature associated with the spatial metric jij. 
The projected derivative operator reduces to the spatial covariant derivative Di associated 
with 7jj, and the extrinsic curvature takes the usual expression in adapted coordinates, 

K^, = ^ {^a^J - D.N, - D,N,) . (3.3) 

The acceleration gets now the 3-dimensional expression 



a' = D'\nN, a = v^7~a^ , (3.4) 

where D^ = Y''Dj (notice that a° = 0). 

The 3d theory (13. 2 p is equivalent to the previous id theory (12.41) . as is explicitly checked 
in Appendix El at the level of the field equations. However the spirit is different, since the 
r-field disappears from (13.20 and is now hidden in the time coordinate. For instance, the 
acceleration's norm a was a r-dependent quantity in (12.40 . but it depends now only on the 
lapse A^ which is a geometrical degree of freedom. 

Hence one may interpret the new a-dependent term in the 3+1 Lagrangian (13. 2p as a 
Lorentz invariance violating term. While the 4d theory (12. 4p is generally covariant, i.e. 
invariant under the full group of 4d diffeomorphisms 

x^ = x^lx") , (3.5) 



[with r being a scalar field, i.e. t{x) = t{x)], the 3d formulation (13. 2p is only invariant under 
the subgroup of diffeomorphisms leaving invariant the preferred time foliation, namely 

i=i{t), (3.6a) 

x' = x'{x^,t). (3.6b) 

Under this subgroup of diffeomorphisms the lapse and shift transform as 

N=%N, (3.7a) 

at 

dt dx^ ,, dx^ dx^ /„ „, X 

while the spatial metric transforms like a spatial tensor. Note that (I3.7ap corresponds in the 
Ad formalism (such that t = r) to a reparametrization of the time foliation by f{T) = t{t). 

The broken diffeomorphism invariance (13. 6 p has recently been used in the context of 
quantum gravity to build a power-counting renormalizable quantum gravity at high energy 



35|-|37] . As is known, the full diffeomorphism invariance can be restored in the Ad formulation 
at the price of an additional structure. This is the role of the Khronon field r, which appears 
to be the Stiickelberg field |42:] associated with the broken diffeomorphism symmetryo 

To write the field equations in the most convenient way we use the following notations 
as a shortcut for describing the matter part of the field equations: 

^ '^^™^ (3.8a) 



^" - ihfk ■ ''■'-' 

These quantities reduce to the usual notions of the energy density, current density and 
spatial stresses in the case of a perfect fluid. However they apply to any kind of matter fluid 
or field. The variation of the Lagrangian (13. 2p with respect to the lapse N yields a modified 
Hamiltonian constraint equation, namely 

n + K^~ KijK'^ - 2/ + Axa^ + 4A(xa') = 167re , (3.9) 

where we recall that both /(a) and x{^) = /'('2)/(2o) are functions of the norm of the 
acceleration given in (13. 4p . The variation with respect to the shift iVj yields the momentum 
constraint equation, i.e. 

Dj {K'^ - Y^K) = -Snr . (3.10) 

This equation is unchanged with respect to general relativity, since the term we have added 
in the Lagrangian does not depend on the shift. Finally, varying with respect to the spatial 
metric 7jj gives the evolution equation 

g^^ + 1a {k^' - Y'K) + ^ A [iv(*(ir^> - 7^>ir)] 



4 



A well-known example of Khronon field is the one introduced in the Cartan theory of gravity which is the 



covariantization of Newton's theory [43|-|4 



- ^ {D'D^N - -i'WkD^N) - Ixa'a^ + j^'^ = SnV^ , (3.11) 

where Q^^ = TV'^ — ^Y^TZ is the 3-diniensional Einstein tensor, and the convenient notation 
Df = dt — N^Dk is used. We gave here the full equation for completeness, but the only- 
difference with general relativity lies in the two last acceleration-dependent terms with /(a) 
and x(a). The trace of (13.111) gives 

7^ + ^DtK + ?,KijK'^ + K^- ^DiD'N - 6/ + 4xa' = -167r (t + j^NiX^ , (3.12) 

where T = 'jijT^-' , and we also used the momentum constraint equation (I3.10p . By elimi- 
nating TZ using the Hamiltonian constraint (13. 9p . one obtains an important equation, which 
plays the role of a modified Poisson-type equation for the gravitational field, 

A [(1 + X)a'] +/ + «'- l^DtK - K'^Ki, = ir^ (e + ^N,J' + t\ . (3.13) 

We check in Appendix |A] the compatibility of the field equations in the 3d and 4d rep- 
resentations, that are obtained by varying with respect to different degrees of freedom, i. e. 
g^u, T and \l/ in the Ad formulation, and A^, Ni, 'jij and \1/ in the 3d formulation. Starting 
with the 4d Einstein equation (12. 5p and applying the projection procedure which allows to 
get the 3+1 equations with the help of the Gauss-Codazzi relations, we get the same 3+1 
equations as obtained by varying directly the action expressed in adapted cordinates, namely 
(I3.9p - (l3.1ip . In other words, starting with the action in 3+1 coordinates and varying the 
non-covariant degrees of freedom, or varying the 4d covariant action and then projecting the 
equations obtained in adapted 3+1 coordinates, yield the same field equations. 

A related question is the fate of the Khronon field equation (12. 9p when switching to 
adapted 3+1 coordinates. Since r becomes the time coordinate, we no longer vary the 
action with respect to it, and the three equations (I3.9p - (l3.1ip do not manifestly contain the 
r- field. In adapted 3+1 coordinates the equation (12.91) becomeqfl 



Dt 



D.,{xa') + Xa'-^ 



+ NK(^Di{xa') + xa') - Nxa'a^K.j = , (3.14) 



where Dt = dt — N^Dj.. We substitute the term Di{xa^) by using the Hamiltonian constraint 
(13. 9 p and perform a series of transformations with the help of the other equations (I3.10p - 
(13. lip . A long calculation, in which one uses the useful relation 



DiD, 



N [K'^ - Y^K) = i A7^ + NK'^n^j , (3.15) 



then shows that (I3.14p is equivalent to the continuity equation for the matter field, i. e. the 
evolution equation of the matter energy density e defined in (I3.8ap . which is given by 

Dts + NKe + Di{NJ') + fD,N + 2KijN'P + NV^ Kij = . (3.16) 



Note that in 3+1 adapted coordinates, we have Vp(xa^) — Di{xa^) + xa 



2 



This equation is the 0-th component of the conservation law for the matter field, i.e. 
y^j^Oi/ _ g_ Furthermore, another long computation shows that, given the 3d contracted 
Bianchi identity DjQ^^ = 0, and the equations fl3.9p - fl3.lip . the Euler equation of the matter 
field is satisfied, i.e. the evolution equation of the current J* given by 

Dtf + NKX + Dj [NV^ + N'J^ + N^ X] + eD'N + JWjN' + 2NXK} = . (3.17) 

This corresponds to the spatial component of the matter conservation Vi/T**^ = 0, once a 
term is substituted using the continuity equation f l3.16p . These checks are the transposition 
in the 3+1 formalism of the fact that the r-equation (12. 9 p is implied by the Einstein field 
equation (12. 5 p (see |38| and Section HTj) . 



IV. NON-RELATIVISTIC LIMIT 

We now investigate the non-relativist ic limit of the model, which will determine the 
required behaviour of the function /(a) in order to recover a MOND-like modification of the 
Poisson equation in the relevant weak field regime. In this Section we restore the G and c 
factors; in the field equations (13.111) for instance, the right-hand side must be multiphed by 
G/c^ so that the stress-energy tensor has the dimension of an energy density. 

We consider a matter system isolated from the rest of the Universe and at rest with 
respect to the preferred frame, in the sense that the fiducial observers have a coordinate 
velocity w^ = n^/n^ = —N^ which is of small post-Newtonian order. We make the usual 
post-Newtonian ansatz on the metric components generated by the system: 

iV = 1 + 4 + 0(4) , (4.1a) 

N, = 0{3), (4.1b) 

7., = '^., (l - f^) + C? (4) , (4.1c) 

where 0{n) denotes a small remainder term of order O (c~"). We introduce in (14. ip two a 
priori different Newtonian-like potentials and ip for the lapse and the spatial metric, but 
we shall check that they are actually identical like in general relativity. 

The acceleration (13. 4p of the fiducial observers reduces to the Newtonian acceleration of 
the 0-potential, 

a=lv0 + O(4) . (4.2) 

We use bold-face notation to represent ordinary three-dimensional Euclidean vectors, e.g. 
a = (a*). With our convention, a has the dimension of the inverse of a length. From (13. 3p 
we see that the extrinsic curvative is a small post-Newtonian quantity, of order Kij = 0{3). 
Finally we assume that the matter source is non-relativistic, so that its energy-momentum 
tensor is dominated by the Newtonian rest mass density denoted p: 

e = pc^ + O{0) , (4.3) 

with comparatively negligible current and stress densities, 

f = pv' + O (2) , (4.4a) 



T'^ = pvV + C (2) . (4.4b) 

Then the equations (13.161) and (13.171) reduce to the standard Newtonian continuity and Euler 
equations (with 0-potential) . 

We first show that the two Newtonian potentials and ip are equaL This follows from 
the trace of the evolution equation given by (13.121) . The function /(a) scales with a^ where 
the dimensionful acceleration is of order a = 0{2) from (14. 2p . This means that the post- 
Newtonian order of this function is /(a) = 0{4). Note that since /(a) will contain a 
cosmological constant A (see e.g. (14.81) below), to be consistent the post-Newtonian order 
of the cosmological constant must also be A = 0{4). Inserting into (13.121) the previous 
estimate /(a) = C(4) together with Kij = C(3), and using the fact that at the dominant 
order we have TZ = -^Aip + 0(4), we readily obtain A{ip — 0) = C(2). For a regular isolated 
matter source the only solution is 

ij = ^ + 0{2) . (4.5) 

The equality of the two potentials is very important for the viability of the theory as an 
alternative to dark matter (see J46[ for a review). It implies that the light deflection and the 
gravitational lensing will be given by the same formula as in general relativity, but with a 
single potential (j) = ip obeying the modifled Poisson equation (given by (14. 6 p below). Unlike 
in Te VeS |l5l . |l6| , there is no need in the present theory to invoke a disformal coupling of 



the photon fleld to gravity in order to get the light deflection viable. The property (14. 5 p has 
been proved in (l9| and is a generic feature of non-canonical Einstein-aether theories. 

We now turn to the equation satisfled by the Newtonian potential 0. Clearly, from the 
form of the metric (14. ip . the ordinary baryonic matter feels the Newtonian gravitational 
fleld g = — V0. From (14. 2 p we have a = —g/c^ at the leading order. In the non-relativistic 
limit, and for the non-relativistic matter source satisfying (I4.3p -( l47^ . the equation (I3.13p 
becomes 

V- (l + x)V0 =47rGp + C(2) , (4.6) 

which exactly reproduces the MOND equation (II. ip for the gravitational fleld g = —'V(j). 
The MOND function is simply given hj n = 1 + x, and x can be interpreted as a "susceptibil- 
ity" coefficient for the modifled Poisson equation. Using the fact that a = g/c^ at the leading 
order, the latter correspondence between fj,{g/ao) and the function xi^i) = /'('3')/(2a) will 
tell us what are the constraints to be imposed on the function /(a) in the initial Lagrangian. 
To recover the Newtonian regime when a ^ oq {i.e. formally when a — )■ oo), this function 
should tend toward a constant. More precisely, we shall recover general relativity with a 
cosmological constant which is the limiting value of /(a), say 

/(a) = Aoo when a -> oo . (4.7) 

On the other hand, to recover MOND asymptotically, we require that when a — )■ (a has 
the dimension of [length]"^), 

2 2c2 3 „ /cV 



f{a) = Ao-a' + —a' + 0{^\ , (4.8) 

3ao V «o / 

where Aq is constant, and we recall that the MOND scale is measured to the value oq — 
1.2 X 10~^° m/s^ y^. We see that the term oc a^ in the expansion (14. 8 p is similar to one of the 
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terms proposed in the "healthy" extension of Hofava gravity |36l. l37l|. Using the notation of 
36l . |37| the coefficient in front of this term in the action is a = 2. Beware, however, that 



this term here is only the first term in a weak-field expansion. 

A priori, the cosmological constant Aqo appearing in the general relativistic limit (14. 7p 
could be different from the cosmological constant Aq in the MOND regime (14. 81) . Notice 
that in an homogeneous and isotropic Universe, making the r-field to coincide with cosmic 
time, a vanishes identically and the equivalent energy-momentum tensor f l2.10p reduces to 
a cosmological constant /(O) = Aq. However, the precise cosmological implications of the 
model would need a study at the level of perturbations, which is left for future work (see 



23H26| for cosmological studies in modified gravity theories). For simplicity, we provisionally 
assume here that Aqo and Aq are equal and we denote their common value by A. 

Note that in the MOND regime at least, the post-Newtonian order of magnitude A = 0{A) 
we have assumed for the cosmological constant is natural. Indeed, because of the numerical 
coincidence between oq and the acceleration scale oa ~ VA. associated with the cosmological 



constant 47|, A scales like clq/c^ so that we indeed have A = 0(4). In this picture we are 
considering that Oq represents a new fundamental constant independent from c and G, 
presumably coming from some new (unknown) fundamental physics. Therefore, factorizing 
out A in (14.81) would yield a Taylor expansion in terms of the small dimensionless parameter 
c^a/ao, with coefficients expected to be of the order of one. 

The coincidence between oq and a^ can be made completely explicit in the following 
example. A function fulfilling all the above requirements is the "Planckian" distribution 



a' 



where the "temperature" oa (in fact kT/^ = hca\) is that of the horizon of the de Sitter 
Universe associated with the positive cosmological constant, namely [48| 



a. ^ yi , (4.10) 

By comparing the expansion when a — )■ of this function with (14.81) we deduce that qq would 
be related to the cosmological constant A by 

ao = -^— . (4.11) 

The numerically value of ao would be extremely close to the actual one: indeed we have 
c^Oa — 0.9 X lO^^'^m/s^ so that Oq — 1.2 x 10^^°m/s^, in very good agreement with obser- 
vations. However, the function (14. 9 p is only an example without physical justification. 

Note that we can also retrieve the deep MOND regime directly at the level of the La- 
grangian (13. 2p . Indeed we have, up to a total divergence, ^J^N1Z = —-^\'Vip\'^ + C'(6), 
in terms of the potential ip parametrizing the spatial metric (I4.1cp . But because the 
two potentials ip and are equal, we see that this term will be cancelled by the term 
2a^ = ^|V0p -|- 0{6) coming from the expansion of /(a) as given by (14. 8p . Therefore 
we find that the Newtonian behaviour is suppressed, and it remains the next term in the 
expansion which is proportional to a^ and immediately gives the MOND Lagrangian 

-Cmond = -T7r4^|V0p - p0 , (4.12) 
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whose variation with respect to reproduces the deep-MOND equation 

\V6\_ " 



ao 



AnGp . (4.13) 



The cancellation of the Newtonian term in the Lagrangian by a "counter-term" oc a^ and 
its replacement by the next-order cubic term oc a^ is similar to what happens in other 
approaches, see e.g. |l9|] and equations (67)-(69) in 32 • 



Thus we have recovered the MOND equation in the appropriate regime. Beware however 
that the analysis presented in this Section relies on the fact that we choose coordinates 
adapted to the foliation. If we come back to a 4d formulation with the r-field being a 
dynamical field, we still may make the hypothesis (14.11) for the metric generated by a system 
at rest in a given frame, but the expression (14. 2 p becomes 



^ = — 

c^ 



1^ + ^^ 



+ C(4) , (4.14) 



with w = n/n^ being the coordinate velocity of fiducial observers in that frame and diu/dt 
being their coordinate acceleration. The fiducial observers will follow a non-geodesic motion 
governed by the r-equation, thus a 7^ in general, but the first term in (I4.14p will appear 
as a preferred- frame effect. Note that this effect depends on the acceleration and occurs 
already at Newtonian order. In the strong field regime where |V0| ^ ao (Solar system, 
binary pulsars, etc.), the presence of this term should not change the fact that a ^ cto/c^, 
so the modification of the dynamics with respect to general relativity should be suppressed. 
In the weak-field regime however, for instance in the outskirts of galaxies, this term could 
lead to a sizeable effect. We would need to investigate the complicated dynamics driven 
by the r-equation to get a precise answer. However if we assume that the preferred frame 
essentially coincides with the frame of the large scale structure or the cosmic microwave 
background, this coordinate acceleration of fiducial observers in a galaxy-centered frame 
becomes (minus) the acceleration of the galaxy with respect to the CMB frame, which is 
typically of the order of magnitude of ao/100 (see e.g. (49|). In this case the term dto/dt in 
(I4.14P is not expected to significantly impact the fit to MOND phenomenology. 



V. CONCLUSION 

We have proposed a relativistic modified gravity model, based on a preferred 3-1-1 space- 
time foliation, reproducing the phenomenology of MOND j^ in the weak- field limit. The 
modification with respect to general relativity consists in adding to its ordinary Lagrangian 
a function of the norm of the acceleration of the congruence associated with this foliation. 
We investigated two different points of view on this theory: In the first, we introduced a 
scalar field called the Khronon defining the foliation by its constant-value hypersurfaces, 
keeping a full 4-dimensional covariant formalism, while in the second, we wrote the theory 
in a 3-1-1 fashion, in a frame where the time coordinate coincides with the Khronon field and 
where the Lagrangian is no longer manifestly Lorentz invariant. The spirit of our approach 
is similar as that of recent attempts at building a consistent quantum theory of gravity 



35H37J. but here the Lorentz invariance violation occurs at low energy, and our motivation 



is purely classical. We showed the equivalence between the 2>d and M formulations of the 
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model, and gave the requirements on the function initially introduced in the Lagrangian to 
recover MOND in the non-relativistic approximation. 

We leave to future work some important questions. First, it is easy to see that in an 
homogeneous and isotropic background, our model simply reduces to the addition of a 
cosmological constant. However the cosmological implications of this model in a perturbative 
regime have not been studied yet. We have in mind the viability of the model when faced 
to cosmological observations, notably the anisotropics of the cosmic microwave background 
and the structure formation (see |23l426l| for general investigations of the cosmology in TeVeS 



and Einstein- ffither theories). 

A second question is that of the viability regarding Solar system tests and binary pulsars 
data. Although it should recover general relativity plus a cosmological constant in the strong 
field regime appropriate to these systems, see fl4.7p rl our model includes a preferred frame 
effect which we would like to quantify precisely, using various hypothesis for the motion 
of the Solar system with respect to the preferred frame. (See |27|, |28| for the computation 



of the preferred frame post-Newtonian parameters in TeVeS and canonical Einstein- a3ther.) 
More importantly, the preferred frame effect should also be computed in the MOND regime, 
using various forms for the MOND function, to see how it affects the usual MOND fit of the 
fiat rotation curves of galaxies; see (I4.14p and the discussion right after. 

Another question is related to the dynamics of the Khronon field itself [obeying the 
equation (12. 9p ]. Indeed in this model the foliation is dynamical, and the evolution of the 
Khronon field should be compatible with the preservation of its smoothness and regularity 
properties allowing a space-time foliation to be built on it. 
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Appendix A: 3-|-l decomposition of the modified Einstein equations 

In this Appendix, we check the equivalence of the 4d and 3d formulations of the theory 
discussed in Sections Hi] and lllli The material used here is well-known, and the derivation 
will basically proceed along the same steps as in standard general relativity (see e.g. |50|), 
but we choose to present it for completeness. 

The projected covariant derivative operator D^ is defined in the usual way by e.g. D^A'^ = 
liflJ'^ p-^ for any vector belonging to the spatial hypersurface, i.e. riuA" = 0. The spatial 
Riemann tensor associated with the derivative D^ is denoted TZ^^^p^. The extrinsic curvature 
tensor K^y = "^^^V pff is symmetric on account of the Frobenius theorem n\^pVunp] = 
satisfied by the hypersurface-orthogonal vector n^. The first Gauss-Codazzi relation 

T^tMupa = Ip^'lflp^l^Rcp-yS + Kp^Kyp - KppKy^ , (Al) 

is contracted with 'jt^P'-f'^'^ which yields 

n = R + 2n^'7fR^, + K^^K^' - K^ . (A2) 



^ For instance, the MOND transition radius for the Sun is ^GAfe/oo ~ 7100AU. 
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Inserting into (1A2|) our modified version of tlie Einstein equations fl2.5p - fl2.6p . we obtain tlie 
modified Hamiltonian constraint as 

7^ + i^2 _ ^^^j^^u _ 2 J + 4v^ (;^a^) = levr n^n^T^, . (A3) 

Next, the second Gauss-Codazzi relation, 

D^ [K^"" - Y^K] = Y^nPRup , (A4) 

when projecting the modified Einstein equation along '-j^'^n'', leads to 

D^ [K^"" - Y'^K] = SnY^n^T^p . (A5) 

The /(a)-dependent terms disappear in this second constraint equation, as expected since 
it corresponds to the variation of the action with respect to the shift, and our modification 
does not concern the shift. Using the fact that 

n'n'^R^pu. = K^pK/ + ^D^D.N - ^C,K^, , (A6) 

one contraction of the first Gauss-Codazzi relation yields 

T^,u = l/lfRp. - ^^iK,u + ^D,D,N + 2K,pK/ - KK,, . (A7) 

Here Ci denotes the Lie derivative with respect to the vector £^ = Nn^. Projecting the 
modified Einstein equation on ^J'^J'-, we find: 

+ KK^, - 2K^pK/ + 7^, i^K^^Kp, - \k^) + h,, - 2xa,a, = 87r7/7/T,. . (A8) 

Finally, making use of the correspondence [compare also with fl3.8p 

e = Ar2T°° , (A9a) 

/ = Ar(r°^ + A^^T°°) , (A9b) 

7-ii = T'^ _ ArW^T°° , (A9c) 

one can check that, when using the adapted 3+1 coordinates, the equations (IA3p . (lASp and 
(1A8I) are the same as the 3+1 equations fl3.9l) - fl3.1ip obtained by direct variation of the 3+1 
Lagrangian (13. 2p . 

Appendix B: Stress-energy tensor of the equivalent matter field 

As we discussed in Section [TTl we can consider the r-field, in the 4d formulation, as an 
additional matter field. Its stress-energy-momentum tensor is then 



rpp.u 

" ~ 87r 



f{a)g'"' + 2n^'n''Vp{x{a)a'') -2x{a)a^'a'' . (Bl) 
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If we choose n^ as the natural 4-velocity associated with this matter fluid, we can perform 
the canonical decomposition of the energy-momentum tensor according to 



Tr = {er + Pr) n^rf + PrQ^" + 2Q[^7f^ + S^ 



^^u 



(B2) 



where Sr and Pr denote the energy density and pressure of the fluid, where Q^^ is the heat 
flow orthogonal to the 4-velocity, i.e. n^Q^ = 0, and where T,!^'^ is the symmetric anisotropic 
stress-tensor which is transverse to the velocity and traceless, i.e. n^T,^'^ = and g^yT,';^^ = 0. 
We then obtain 



1 



6T[ 



Pr 

Q'r 



1 

Stt 
0, 

x(q) 

47r 



2 



-/ + la\ 



a^'a" 



-7'^^a2 



(B3a) 

(B3b) 
(B3c) 

(B3d) 




In terms of those definitions, the 3+1 version of the conservation equation '^ JT!^^ 
reduces to a scalar equation, taking the form of the "continuity" equation 

DtEr + iV (e, + Pr) K + NYi'^Kij = . (B4) 

As we mentioned in Sections [III and IIIH the conservation equation V^T^'^ = has exactly 
the same content as the Khronon field equation f l2.9p . and is equivalent, via the modified 
Einstein equations f l2.5p . to the conservation of matter V JT^^ = 0. Thus, this equivalent 
fluid has the property that its Euler equations reduce to its continuity equation, i.e. the 
four components of VyT!^" = reduce to only one equation (as expected for a scalar field). 
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